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I Intioduction 

This report is the final report for the Cooperative Agreement NCC-2- 
1000. The tasks outlined in the various proposals are listed below 
with a brief comment as to the research performed. The publications 
are attached. 


II. NCC 2 1000 Tasks and Accomplishments 

1. Development of an aerosol chemistry model (Tasks 1-5, Dr’s 
Tabazedeh and Lin) 

This was the main focus of this Cooperative Agreement and the 
results are contained in two published papers attached) 

2 Utilization of satellite measurements of trace gases along with 

analysis of temperatures and dynamic conditions to understand ice 
cloud formation, dehydration and sedimentation in the winter polar 
regions. (Task 6 Dr Stone) 

This was a specific task for Dr. Stone and resulted in a paper has been 
submitted for publication (attached) 

3. Comparison of the HALOE and SAGE 11 time dependencies of the 
Pinatubo aerosol decay (Task 7, Dr. Bergstrom) 

This was a specific task for Dr. Bergstrom and resulted in a paper that 
is being prepared for publication. 
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Abstract, Recent laboratory data on ice freezing from 
aqueous H 2 SO 4 solutions are used to update our homogeneous 
ice freezing nucleation algorithm. The effect of particle size on 
the ice freezing process is demonstrated by calculating and 
comparing ice freezing curves for different size particles 
ranging from 0.2 to 4 micron in radius. Using the ice 
nucleation model we show that liquid water saturation is 
required above -44°C to activate submicron H 2 SO 4 particles 
into cloud droplets. A thermodynamic model is used to show 
that the available laboratory data on ice freezing from ternary 
H 2 SO 4 /HNO 3 /H 2 O solutions are insufficient to adequately 
address the effect of HNO3 uptake on the homogeneous freezing 
of aqueous H 2 SO 4 particles into ice. 

Introduction 

Ice clouds play important roles in Earth's climate 
[Ramanathan et al, 1983] and chemistry Solomon [1990]. 
One mechanism that has been suggested to initiate ice cloud 
nucleation in the atmosphere is homogeneous freezing of 
aqueous H 2 SO 4 solution droplets and Dodd, 1989; 

DeMott et aL, Jensen et al., 1994; Heymsfield and 

Miloshevich, 1995; Demott et a/., 1997]. However, the ice 
freezing process through this mechanism is complicated by 
the fact that a significant fraction of gas phase HNO3 will 
partition into sulfate aerosofs near the ice freezing point of the 
solution [Tabazadeh et ai, 1994; Carslaw et ai, 1994]. 

The ice freezing curve in aqueous H 2 SO 4 droplets was 
recently measured by Koop et al. [1998]. The ice freezing 
temperatures determined by Bertram et al. [1996] for aqueous 
H 2 SO 4 droplets, which we previously used to parameterize our 
ice nucleation model [Tabazadeh et ai, 1997 a, b], appear to be 
systematically loo high [Koop et al, 1998]. In addition, 
differential scanning calorimetry has recently been used to 
determine the ice freezing curve in aqueous HNO3 and ternary 
H 2 SO 4 /HNO 3 /H 2 O droplets [Chang et al, 1999]. In this work 
we use recent laboratory measurements on ice freezing from 
binary and ternary solution droplets [Koop et al, 1998; Chang 
et al, 1999], along with new surface tension data on the 
aqueous H 2 SO 4 system [Myhre et al, 1998], to update our 
homogeneous ice freezing nucleation algorithm. 
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Model Update 

Figure 1 shows critical ice saturations expected to nucleate 
ice in aqueous H2SO4, HNO3 and ternary H 2 SOj/HNO,/H 20 
solutions as a function of the ambient water vapor pressure 
based on recent laboratory measurements. The data of Koop et 
al. [1998] and Bertram et al [1996] show very different 
behaviors for the freezing of ice from aqueous H 2 SO 4 solution 
droplets. The main objective of this paper is to update our ice 
nucleation model to agree with more recent measurement.s. 

An interesting feature of Figure 1 is that the ice freezing 
curves for both aqueous H2SO4 and HNO3 solution droplets are 
nearly identical. According to equilibrium model calculations 
[Tabazadeh et al, 1994; Carslaw et al, 1994], the ionic 
strengths of H 2 SO 4 /H 2 O and HNO 3 /H 2 O at the same relative 


- H^SO^ (Koop et al.) — - H^SO^ {Bertram et al ) 



H^O Vapor Pressure (mb) 

Figure 1. The variation of critical ice saturation as a function 
of water partial pressure over the solution. The solid and dotted 
curves are calculated directly from the tabulated relations given 
in Koop et al [1998] and Chang et al [1999]. The numbers on 
the dotted lines give the H 2 SO 4 content of a ternary 
HNO 3 /H 2 SO 4 /H 2 O aerosol in w’eighl percent. The solid dotted 
curve is based on the data of Bertram et al [1996] converted to 
an appropriate form using the method described in Tabazadeh 
et al [1997a]. Also shown are the results of our ice nucleation 
model for a 4 micron aqueous H 2 SO 4 particle. 
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humidity arc similar in magnitude. Since the equilibrium 
freezing point depression is often directly propi^rtional to the 
ionic strength, the fact that the tuo binar\ electrolytes ha\e 
similar ice freezing kinetics is at least partialK rationalized. 

Also shown in Figure I are the ice Iree/ing curves for 
H^SOj/HNOi/H T) ternary .solutions, which contain 1 and 5G 
H.SO 4 by weight \Chan^ et al.. 19991. The slight (< ) 

increase in the critical ice saturation in a ternary solution, as 
compared to the binary systems, is rationalized by the increase 
in entropy of the solution as a result of mixing. Below we use a 
thermodynamic model to show that the current laboratcuy data 
on ice freezing trom ternary solututiis are insufficient to 
adequately address the effect of UNO, uptake (by sulfate 
aerosols) on the ice formation process in the atmosphere. 

Denoting (cm') and as the volume and weight percent 
of a H.SO 4 solution droplet, the rate ot ice nucleation (particle 
sec’’ ) at a given temperature is: 


molecules in solution based on laboratory ice freezing diita. 
We tolltnv the s.ime procedure here to derive a new relation lor 
the activation free onergv ol Water molecules in solution based 
on the ice freezing data ot Koof) ci a!. \ 199b]. 

fn deriving an expression h'lr the activation free energy, we 
generated two fittings ol the activation energv tunction for 
particle sizes of 1.55 and 6 ..s micron in radius that bracket the 
labor»Uory size range \Koop t ! a!.. 1998], The final activation 
energy relation reported in Table 1 is obtained bv averaging 
the two expressions above tor different particle size 
assumptions. In addition, we replaced the surface tension 
relation in the model bv a new expression based on recent 
laboratory measurements by M\hrt’ ct al. ||99S|. All the 
functions required for the nucleation calculations are described 
and summarized in Table 1. The nucleation model calculates a 
unique solution for the ice freezing temperature of a given size 
droplet as a function ot onl> the ambient water vapor pressure. 
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where C ( particle sec ‘ ) is the preexponeniial factor, 

(ergs) is the Gibbs tree energy for the formation of the ice 
germ. (ergs) is the diffusion activation energy of water 

molecules across the ice/sulfaie solution phase boundary, and 
k is the Boltzmann constant. The preexponential factor and 
the Gibbs free energy are estimated as follows: 

C( T. U-, . V']; ) = 2. 1 X 1 0-'-’ VC ( 2 ) 


=^na, j (3) 

where is the interface energy between the ice/sulfaie 

solution and r„ (cm) is the critical germ radius. For our 
purposes, an ice freezing event occurs when ./ = 1 particle sec b 
W'e previously described [Tahazadih et al.. 1997 a. bj how 
the equations ( 1 ) through (3) can be used to extract a numerical 
relation for the variation of activation free energy of water 


Atmospheric Implications 

Figure 2 shows the efiect ol particle size on ice nucleation 
from an aqueous H^SO, solution droplet. The size-dependence 
calculated by our model is in good agreement with the size- 
dependence measured in the laboratory for ice nucleation in 
pure supercooled water droplets (.see Figure 2 in Fruppacher 
[1995]). As expected, due to smaller volumes ( \ ^ term in 
equation 1 ), higher ice saturations are needed to form ice in 
submicron atmospheric particles as compared to micron-sized 
aerosols studied in the laboratory \Koop et al.. 1998]. Recent 
measurements [Chen ei al.. 2000] support our calculations 
that submicron sulfate particles freeze at lower temperatures as 
compared to mi cron- si zed particles. 

Laboratory freezing data on submicron water droplets 
indicate that such droplets can often remain supercooled to 
temperatures near and sometimes below -40 C \Ha^en et ai. 
1981: Fruppacher. 1995] Thus the fact that our ice nucleation 
model calculates a liquid water activation temperature below 
-40 'C for submicron aqueous particles (Figure 2) is not 

surprising since the efiect of solute is to further depress the 
freezing point of ice in solution beyond that of a pure 
.supercooled water droplet. Overall, a nucleation model offers 


Table 1. I ce Nucleation Functions 


Latent heat of melting in erg mol ' { T is in K): ( T) =10 |6005 2356 ■+ 1 8,27 1 9( T - 273. 1 5} - () ()6.^.*'-^ T ~ 273 1 5 e J 

Ice densiiy (7' is in K) ; p„,(Tl = 0.916- 8-75 x lO'S? - 273, 15) - 1 .667 v IO‘S/ -27.3.15)- 

Sulfate soluiion/air interface energy based on the data of Myhre et ai [ I998j in dyn cm ' (where u is the H.SOj weight percent of the solution): 
= 85.755071 14 + 9..5419663I8.V10’- u - 1 , 103647657.vl 0'‘ >i ' + 7.485866933 a lO'Sv' - 1 .91 2224 1 54,vl0" u '’ + I 7.'6789787.> I O'Ui '' 


^ ,\ul/ oir 


= S2.01 197792 + 5.3I2O72O92.VI0"' i. - 1 ,050692 1 23.vlO'' n - +5.4 1,32606 1 7.U0”’u -' - l -1455"3S27 vl0"' « ■■ + 8,969257061,»10 ' vc' 
cr-^l . = 77 406S2664-6.963l23274,vl0‘Si -9.682499074.vl0'Si- + 8,87979S80.vI0'^ h"’ - 2.3K4^6951 6.^U^^^ ■“ + 2.095358048.»10'Ui- 


For 180 < 7- < 220 K. + (5.5 - 0.0257 )(ff;,^);„„ ■ 

For 220 < 7- < 260 K. +(6.5 - 0.0257- 




Diffusion activ ation energy lO'^ in ergs derived from the Koop et ol. [ 1998] ice nucleation data ( / i> in K); 

ForT < 220 K. = -17459.5l6l83 + 458,4582755l7'-4 849283I3I7'7'- + 0.026(K)3658878r ' - ‘ ;99l57776yxlO' 

+8,9049094618.vl0'’' 7'" - 2.49,32257419,1 1 0‘" 

For T> 220 K. A/-;„, = 104525.9.3058 - 1 103.7644651/ + 1 ,0703327027-- +0,OI7386254322r' - 1 55"685426S.\ lO”' ■/' 
-3.26619I2497.V10-’'-/--' + 6.4679544,59.vl0-"'/-'' 

Equilibrium sulfuric acid weight percent composition and water saturation vapor pressure: See Tahaziuiih et al [ 1997a] 
Sulfuric acid solution densiiy: See Myhre et al. [ 1998] 
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H^O Vapor Pressure (mb) 

Figure 2. The variation of the crincal ice saturation (and 
supercooling) as a function of water partial pressure and 
particle size in radius calculated u^lng the ice nucleation 
model. The critical supercooling is defined as the difference 
between the critical ice freezing temperature and the ice frost 
point of the solution. Liquid water acti\ation temperature^ 
marked on the plot show transition temperatures above which 
particles activate into liquid water droplets prior to freezing 
into ice. The variable X in the above equation is In (PH 2 O). 
where PH.O is the ambient water vap^r pressure in mb. 


several advantages compared to parameteri/atii>ns for treating 
the ice formation process in the atiruvsphcre. Simple 
parameteri/ations tend to provide information that is 
exclusive to one size panicle, whereas in the atmosphere ice 
nucleation in principle occurs over an entire submicrun 
aerosol particle size distributu>n. 

laquid water activation temperatures marked in Figure 2 tfor 
three different particle si/es) shnw iransiiion tempeiaiurcs 
above which sulfate particles will first activate into liquid 
water droplets prior to freezing into ice, To accurately predict 
this transition point in the atmosphere is important since 
cirrus ice ch>iKN that liisi go through a liquid water activation 
phase are often composed ol a large population of small ice 
crvstals t'( al., 199Sj. Radiauvu properties of emus 

clouds are strx'ngly aflected by the numhei l^l lee particles 
luiclcaied \Jcnwn anJ /'(xni, 'Fhe results of our 

calculations siiovs that suhnncnHi II SO; panicles in the 
aimo>{)hcrc can remain sufiei cooled lo temperatures near -44 
C 'I'herefore u is likely that above -44 (' cirrus icc clouds 
olten form bv rtrsi going through a liquid water activation 
pluisc provided that they arc devoid of heterogeneous solid 
impurities Demon ei o/. . 1997). 


The composition of an aqueous H^SOj solution droplet in the 
atmosphere is altered by the uptake ot HNOv Li Figure 3 the 
change in the ternary solution composition as a tunction ot 
temperature is shown for three selected water vapor pressure 
profiles, which bracket typical conditions in the lower 
stratosphere and upper troposphere. The shaded area shows the 
range of H^SO^ compositions of ternary systems tor which the 
ice freezing properties of solutions have been determined in 
the laboratory [Chan^ ei al., 1999], 

Ternary solution droplets containing betw'cen 1 to 5 w't 9c 
H>S 04 freeze at ice saturations near that ot the binary systems 
[Chang tn ai, 1999, see Figure 1). In the background 
stratosphere the H.SOj weight percent in solution near the ice 
nucleation point is located wiiliin the gray area for which the 
ice freezing properties of ternary solutions have been studied 
in the laboratory. Thus the effect of HNO;, uptake in depressing 
the ice freezing temperature is negligible (Figure 1) it the 
stratosphere is in a background state (Figure 3). However, tor a 
volcanically perturbed state, where ILSO^ mixing ratios are 
higher, the H.SOi weight percent of a ternary solution near the 
ice freezing point is close to 15 9- (Figure 3), Because this 
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Figure 3. 4’he variation of lernaiy HjSOj/HNCO ;/iLO aerosol 
composition with lempeialure for diMereni water vapor 
pressure regimes based on the thernuKh namic relations given 
in Luo et a!. \ I995|, The 0 I iind O.Ol mb curves are calculated 
h)r a constant tLSCO. atmospheric mixing ratio o{' 1 00 pp[ and 
variable H.\(); atmospiienc mixing ratios ot 0. 1 i backgrrmnd) 
and 2 pph (perturbed hv convective pol Union) at 200 mb. The 
0.0003 mb curves 4 ire Cvilculated lor a constant UNO;, 
atmospheric mixing r.itu* of lO pph and vanahle M-SO 4 
mi.ving ratios ot (,).3 i background > and 20 pph f perturbed by 
volcanic emplumsi vtt Uf) mb. TTie calculaiioii'' were [leilorined 
down to about -3 2 K t'clow ihe eqmhbruiin ice condensation 
poini (0 the soUmon. the slualed area shows the range ol 
ILSO; coin[>osinons in lemarv aerosol solutions tor which the 
ICC tree/ing properties ot ternary systems have been measured 
m the lahoraloiv |fV/o7?c a o/.. lOVdj. The anows m the plot 
ivler to tl^SOi composition 111 (eriurv solutions tor which 
laboratory ice trec/ing data are lacking. Lahoiaiory freezing 
data \\)T lernarv sttlutions are cuneiulv uiuivaiiahle lor vmicanic 
and background conditions m the stratosphere and upper 
troposphere, respecti veiv . 
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ternary composition regime has not been studied in the 
laboratory, it is difficult to compare and contrast the effect of 
HNO 3 uptake on the ice freezing proce>> in the stratosphere 
between the two cases. Satellite data indicate that it may be 
more difficult to freeze aerosols in a olcanically perturbed 
environment as compared to a background state [Santee et al.. 
1998). Thus additional laboratory data are needed to better 
understand how volcanic eruption> may affect the ice 
formation process in the lower stratosphere. 

Similarly in the upper troposphere the laboratory data on ice 
freezing form ternary systems are insufficient to adequately 
describe and compare the differences between the background 
and perturbed states. For example, in a polluted upper 
troposphere HNO; mixing ratios reach values as high as 2 
ppbv [Laaksonen et al.. 1997; Schnieuer et al.. 1998], The 
high levels of HNO 3 in the upper iropo>rhcre are linked to dry 
convection that injects polluted boundary layer air directly 
into the upper troposphere. The H;;S 04 weight percent of 
ternary solutions in a polluted upper troposphere are located 
within the gray shaded area (see Figure 3). On the other hand, 
ice freezing curves for a background state (between 8 to 12 ^ 
H 2 SO 4 and 1 to 5 9c HNO 3 by weight >. containing about 100 
pptv of HNO,, have not yet been explored in the laboratory 
(Figure 3). Thus it is difficult to address the question of how 
convection (leading to high levels of HNO,) may affect cirrus 
formation in the upper troposphere when no information is 
available on ice freezing properties for compositions 
representative of the background state (--1 00 pptv of HNO,). 

Summary 

We updated our homogeneous ice freezing nucleation code 
using recent laboratory measurements. Our results indicate that 
ice freezing curves for submicron aerosol particles are 
significantly different than those determiined in the laboratory 
for micron-sized particles. Further we show that in the 
atmosphere, for temperatures abo\e -44 C. submicron aquetius 
H.SO 4 particles will activate into liquid water droplets prior to 
freezing into ice. Finally the available laboratory data on ice 
freezing from H^SOj/HNO^/H-O solutions are insutticicnl to 
address how HNO, uptake by aqueous H^SO^ particles may 
affect the ice formation process in the atmosphere. 
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Abstract. Simple expressions are fitted to the results obtained from ion interaction 
thermodynamic models for calculating HNO 3 and H 2 O vapor pressures over 
the NH 3 /H 2 SO 4 /HNO 3 /H 2 O system at cold temperatures. The vapor pressure 
expressions are incorporated into a mass conserving equilibrium solver for computing 
aerosol compositions in the lower stratosphere and upper troposphere. The 
compositions calculated from the aerosol physical chemistry model (APCM) are 
compared against previous pararneterizations. The APCM compositions are in 
better agreement with the compositions obtained from ion interaction models than 
from other previous fonniilations of the NH 3 /H 2 SO 4 /HNO 3 /H 2 O system. The only 
advantage of the APCM over the ion interaction approach is that the numerical 
scheme used in the model is fast and efficient for incorporation into large-scale 
models. The APCM is iis( 3 d to calculate HNO 3 solubility in ammoniated aerosols 
as a function of HNO 3 . H 2 SO 4 and NH 3 mass loadings in the lower stratosphere 
and upper troposphere. While the uptake of HNO 3 by ammoniated aerosols is 
strongly dependent upon the solution neutrality (or pH), we find that in both 
the lower stratosphere and upper troposphere a significant fraction of HNO 3 will 
exist in aerosol solutions near and below the ice frost point irrespective of solution 
neutrality. 


1. Introduction 

Therinodynainic electrolyte m(;dels are often used for 
calculating properties of inorganic aerosols in the lower 
troposphere [Stelsoji and Seinfeld. 1981; Pilijiis and Se- 
infeld, 1987; Wexler and Seinfeld. 1991; Kirn and Se- 
infeld, 1995; Jacobson et ai, 199G; Jacobson^ 1999b] 
and stratosphere [Tabazadeh et ai. 1994; Carslaw et ai. 
1994, 1995b; Weisenstein et al., 1997]. In the lower tro- 
posphere. aerosol models are often used iti air quality 
studies to iissess the effects of aerosols on health, gas- 
phase partitioning, and visibility. In the stratosplu‘re. 
aerosol models have been used to simulate the forma- 
tion and growth of polar stratospheric clouds, which are 
linked to stratospheric ozone depletion [Solomon. 1999). 
However, thermodynamic treatments in large-scale at- 
mospheric models are in general not suited for calcu- 
lating aerosol compositions in the ui)per troposphere. 
Since upper tropospheric aerosols [;artici[)ate in the nu- 
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cleation and growth of cirrus clouds and may also be 
involved in the scavenging of trace gas species, it is 
important to understand their chemical and physical 
properties [Kdrcher and Soloinon, 1999]. 

Recently, equilibrium aerosol formulations have been 
incorporated into three-dimensional models to simulate 
the radiative impacts of aerosols on climate [Adams et 
ai, 1999; Jacobson, 2000]. In this work, an aerosol phys- 
ical chemistry model (APCM) with an efficient solv- 
ing scheme suitable for incorporation into large-scale 
atmospheric models is developed. The APCM is com- 
pared against various pararneterizations. including re- 
sults from the aerosol inorganics model (AIM2) of Clegg 
et ai [1998a]. 

2. Background on Thermodynamic 
Aerosol Models 

2.1. Ion Interaction Approach 

The ion interaction approach is originally bcused on 
the work of Pitzei' [1991]. The solution behavior in 
Pitzer’s method is determined by a series of single- 
ion (and water) activity equations that are based on 
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Table 1. An Coefficients for the Vapor Pressure of HNOa^ 


r 



^2 

a3 

a4 

as 

ae 

a? 

0.0 

12.6923 

- 0.0325979 

- 0.071719 

0.812331 

2.34482 

0.000438486 

- 0.00428908 

- 0.000709075 

O.l 

12.4691 

- 0.0019955 

- 0.077484 

0.727085 

2.12858 

0.000259723 

- 0.00442994 

- 0.000706633 

0.2 

12 . 2.326 

0 . 026093.3 

- 0.085355 

0.644593 

2.52519 

9 . 32848 e -05 

- 0.00453571 

- 0.000694942 

0..3 

11.9840 

0 . 0.513596 

- 0.094561 

0.567113 

2.63058 

- 5 . 53402 e -05 

- 0.00461276 

- 0.000678339 

0.4 

11.7205 

0.0725417 

- 0.104727 

0.500103 

2.74278 

- 0.000172282 

- 0.00466394 

- 0.000659102 

O..5 

11.4289 

0.0880553 

- 0.116820 

0.450799 

2.86783 

- 0.000241534 

- 0.00468708 

- 0.000632214 

0.6 

11.0837 

0.0983196 

- 0.134271 

0 . 41.5462 

3.02666 

- 0.000267273 

- 0.00467490 

- 0.000579027 

0.7 

10.6697 

0.106923 

- 0.160517 

0.378195 

3.24022 

- 0.000284624 

- 0.00461980 

- 0.000480710 

0.8 

10.2087 

0.1 17295 

- 0.194708 

0.324242 

3.50279 

- 0.000327867 

- 0.00451884 

- 0.000342092 

0.9 

9.73753 

0.130202 

- 0.232997 

0.251054 

3.79040 

- 0.000406493 

- 0.00437409 

- 0.000184569 

1.0 

9.28445 

0.144365 

- 0.271608 

0.165018 

4.07962 

- 0.000510598 

- 0.00418932 

- 2 . 90206 e -05 

1.1 

8 . 87 . 3.35 

0.156600 

- 0.306440 

0.080649 

4.34517 

- 0.000611458 

- 0.00396949 

0.000101871 

1.2 

8.55964 

0.157392 

- 0..327751 

0.039034 

4.52665 

- 0.000618013 

- 0.00372273 

0.000153819 

1..3 

8 . 3.3474 

0.134427 

- 0.331971 

0.092212 

4.60197 

- 0.000411043 

- 0.00345255 

0.000106943 

1.4 

7.95881 

0 . 1 L 0127 

- 0.354342 

0.142888 

4.78919 

- 0.000208202 

- 0.00313268 

0.000158658 

1.5 

7.51749 

0.0956133 

- 0.390173 

0.143658 

5.058 15 

- 0.000120908 

- 0.00275735 

0.000282545 

1.6 

7.14929 

0.0860305 

- 0.421553 

0-118501 

5.29796 

- 0.000104343 

- 0.00234767 

0.000379673 

1.7 

6.88502 

0.0761817 

- 0.441786 

0.092775 

5.46818 

- 0.000107787 

- 0.00192436 

0.000413991 

1.8 

6.68262 

0.0630763 

- 0.453078 

0.080818 

5.58436 

- 0.000101314 

- 0-00149773 

0.000397780 

1.9 

6.43124 

0.04 55338 

- 0.465952 

0.086457 

5.71507 

- 7 . 48198 e -05 

- 0.00106043 

0.000386668 

2.0 

5.72588 

0.0242624 

- 0.524014 

0.096528 

6.13955 

- 5 . 3 1470 e -05 

- 0.00053951 

0.000609824 

OC^ 

11.5869 

0.222271 

- 0.109389 

0.279027 

2.83237 

- 0.00101574 

- 0.00488557 

- 0.000704950 

^ A 

= ao + ai Wi ao Wo 

-f -f- a,4>/VN -f a,5W^ -f agwi W2 + 




^ The oc represents the ternary system of NH3/HNO3/H2O (i.e., H2SO4 — 0). 


thermodynamic properties of mixed solutions. Physi- 
cal parameters for the fundamental activity relations in 
an ion interaction model are fitted to laboratory mea- 
surements of solute-water mixtures of interest. Ther- 
modynamic electrolyte models based on the ion in- 
teraction approach have evoh’ed significantly in recent 
years by Clegg and coworkers [e.g., Clegg and Brim- 
blecombe, 1990, 1995a, 1995b: Carslaw et al, 1995a; 
Clegg et ai. 1998a, 1998b], Since the models of Clegg 
and coworkers, collectively referred to as aerosol in- 
organics model (AIM2), are fitted to laboratory mea- 
surements conducted over a broad temperature range 
(< 200 to 328^K), they are more accurate than tho.se 
based on the common binary activity approach de- 
scribed in section 2.2. However, tlie complex nature 
of the ion interaction ap|)roach makes AIM2 computa- 
tionally impractical for three-dimensional applications. 
Here we use AIM2 to generate solution compositions 
and vapor pressures for the NH3/H2SO4/HNO3/H2O 
system for a wide range of humidities and temperatures. 

2.2. Binary Activity Approach 

The second ap|)roach first ertierged from modeling 
inorganic aerosols in air quality studies [e.g., Stetson 
and SeinfeUL 1981; Bassett and SemfeUU 1983; Saxena 
et ai, 1986]. This method separates out water (via 


either the Gibbs-Duhern equation or the water equa- 
tion and Robinson, 1966]) and solute activity 

coefficients. Lsually. a mixing rule of either Bromley 
[1973] or Kusik and Meissner [1978] is used to estimate 
the mixed solute activity coefficients. The equilibrium 
models based on this approach require only a knowl- 
edge of water and solute activity coefficients at a binary 
level. Thermodynamic models based on the binary ac- 
tivity approach are computationally more efficient but 
less accurate than the ion interaction approach because 
the physics of the activity coefficients in the former 
approach are mainly based on the behavior of binary 
solutions instead of mixed solution properties used in 
the latter formulations. To contrast tw'o different ther- 
modynamic tn^atments, we will update an equilibrium 
model based on the binary activity approach, EQLT- 
SOLV II [Jacobson et a/., 1996: Jacobson, 1999b], and 
compare its predictions against results obtained from 
AIM2. EQLTSOL\' II applies a well-converged numeri- 
cal solver scheme to simultaneously solve a large number 
of equilibrium equations. 

2.3. Vapor Pressure Approach 

A combination of the two approaches outlined in sec- 
tions 2.1 and 2.2 is used in this work to develop a 
fast and accurate parameterization of the NHj/HoSO^/ 
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Table 2 , Bn Coefficients for the Vapor Pressure of HNOa^ 


r 

bo 

bi 

bo 

bs 

b4 

b5 

be 

b7 

0.0 

-7672.46 

67.7280 

32.,3655 

-■259.196 

-235.30 

-0.256884 

0.97336 

0.197376 

0.1 

-7609.77 

50.7068 

32.9809 

-222.190 

-255.43 

-0.131368 

1.07699 

0.201805 

0.2 

-7541.22 

35.7444 

34.4801 

-189.780 

-280.54 

-0.026085 

1.16451 

0.202182 

0.3 

-7467.24 

22.7871 

36.6718 

-162.151 

-309.62 

0.059118 

1.23812 

0.199536 

0.4 

-7387.11 

11.9765 

39.4988 

-140.155 

-342.40 

0.122744 

1.29896 

0.194190 

O.o 

-7297.89 

3.4748 

43.2297 

-124.533 

-,380.56 

0.163341 

1.. 34688 

0.184702 

0.6 

-7194.29 

-3.0516 

48.6478 

-113.773 

-428.87 

0.184446 

1.. 38064 

0.166777 

0.7 

-7074.27 

-8.6193 

56.4190 

-103.509 

-491.39 

0.196255 

1.39907 

0,136740 

0.8 

-6944.62 

-14.0870 

66.1301 

-90.2.35 

-565.46 

0.207.582 

1.40222 

0.096867 

0.9 

-6815.18 

-19. .5849 

76.6808 

-73.739 

-644.27 

0.220387 

1.39134 

0.053245 

1.0 

-6693.25 

-24.6807 

87.0206 

-56.230 

-721.33 

0.231211 

1.36800 

0.011651 

l-l 

-6584.97 

-28.3719 

96.0022 

-42.304 

-789.58 

0.230911 

1.33398 

-0.021.582 

1.2 

-6504.53 

-27.8717 

100.995 

-44.004 

-832.71 

0.192747 

1.29170 

-0.031795 

1.3 

-6449.42 

-19.3426 

100.908 

-77.488 

-843.95 

0.079582 

1.24188 

-0.012861 

1.4 

-6347.58 

-8.4904 

106.038 

-U7.723 

-886.93 

-0.052980 

1.17.584 

-0.022354 

1.5 

-6203.79 

-0.1957 

118.040 

-143.506 

-971.43 

-0.1,56090 

1.08786 

-0.069346 

1.6 

-6055.69 

4.8272 

133.047 

-1.52.440 

-1073.17 

-0.221191 

0.98027 

-0.1.32468 

1.7 

-5927.18 

6.8.545 

148.101 

-146.787 

-1173.72 

-0.249961 

0.85834 

-0.195765 

1.8 

-5829.77 

6.2290 

161.504 

-129.145 

-1262.57 

-0.24.5338 

0.72779 

-0.2,50386 

1.9 

-5768.38 

3.1411 

172..326 

-101.381 

-1.3.34.05 

-0.209351 

0.59425 

-0.291726 

2.0 

-5744.66 

-2.5250 

180.115 

-63.965 

-1,385.49 

-0.141562 

0.46318 

-0.317848 

oc" 

-7298.14 

-51.7642 

35.6176 

-55.070 

-.360.093 

0.376735 

1.41665 

0.252861 

=‘Bn 

= bo + bi 

w 'l -h b ow 

2 T b3 W 1 -j- b4 W' 2 

+ bswf 

+ bew'i w-2 

-h hjW'i. 



^ The oc represents the ternary system of NH3/HNO3/H2O (i.e., H2SO4 — 0 ). 


HNO3/H0O system for atmospheric applications. The 
vapor pressure approach requires generation of large 
sets of HjO and HNO3 vapor pressure arrays over the 
NH3/H0SO4/ HNO3/H2O system for many assumed 
compositions and temperatures. The vapor pressure, 
composition, and temperature fields generated are fitted 
into simple mathematical expressions similar to those 
derived by Luo et al. [ 1995 ] for the H2SO4/HNO3/H >0 
system. The vapor pressure relations are then in- 
serted into a numerical mass- conserving equilibrium 
solver similar to that used in EQITSOLV II [Jacobson 
et a/.. 199 G] for gas-aerosol phas(^ I)artitioning calcula- 
tions. 

3. Model Development 

3 . 1 . Aerosol Physical Chemistry Model 
(APCM) 

Luo et al. [ 1995 ] have shown that etiuilibrium par- 
tial pressures of HNO3 and HoO ovov the H-jSO i/HNO;)/ 
H-jO system roughly follow a Clausius-Clapeyron rela- 
tion of the form InP = A -j- B/T (wheu'c P is pressure. 
T is temperature, and A and B are constants) for a 
fixed solution composition. For the quaternary system 
of NH3/H.)S0.4/HN03/H.)0 we found a similar behav- 


ior for the variation of HNO3 and H2O vapor pressures 
over the solution. 

Assuming that NH3 and H2SO4 reside completely in 
the condensed phase (i.e., both NH3(gj and H2S04(g) 
are negligible), w'e follow the approach of Luo et al. 
[ 1995 ] and define w^ and W9 as the weight percents of 
aminoniated sulfate and nitric acid, respectively, in tiie 
solution as follow's: 

Wi = weight % of (NH4)rH2-rS04(aq). 

wj = weight % of HN03(aq), (I) 

_ NH3 _ >2) 

^-H2S04“hso-^, + so:;-^,^ 

whei'(‘ r in wj is th(^ fixed mole ratio of ammonia to sul- 
furic acid, and (aq) is the aqueous-phase species. The r 
ratio can l>e considered as the degree' of the ammoniated 
solution neutrality and can take on any values in APCM 
(including fractions) between O.U and 2 . 0 . For example, 
if r equals 0.0. 1.0, and 2.0, \\{ represents the weight 
percents of H2SO4, (NH4)HS04. and (NH4)2S04. re- 
spectively, in the solution. Using Cleyg et al. s [ 1998 a] 
model, W'e first generated a series of arrays of HNO3 and 
H2O vapor pressures for a wide range of weight percent 
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Table 3. Cw Coefficients for the Vapor Pressure of H2 0^ 


r 

Co 

Cl 

C 2 

C3 

C 4 

C5 

C 6 

C 7 

0.0 

22.7490 

0.042481’ 

0.0533280 

-0.05674.32 

-0.276555 

-0.000621533 

-0.000311769 

-0.000283120 

0.1 

22.7669 

0.0310031 

0.0493720 

-0.0270045 

-0.265333 

-0.000496930 

-0.000261912 

-0.000257919 

0.2 

22.7838 

0.0223810 

0.0457014 

-0-00o9r'2l 

-0.254223 

-0.000403349 

-0.000224673 

-0.000234421 

0.3 

22.8005 

0.0165207 

0.0423183 

0.0065020 

-0.243407 

-0.000339520 

-0.000198990 

-0.000212594 

0.4 

22.8179 

0.0132215 

0.0393297 

0.0110.504 

-0.233647 

-0.000303073 

-0.000184277 

-0.000192973 

0.5 

22.8367 

0.0119853 

0.0369883 

0.0096303 

-0.226435 

-0.000288721 

-0.000180339 

-0.000177085 

0.6 

22.8556 

0.0119005 

0.0354244 

0.00586.33 

-0.222617 

-0.000287316 

-0.000186497 

-0.000165709 

0.7 

22.8706 

0.012233d 

0.0342757 

0.0024714 

-0.220187 

-0.000291401 

-0.000201103 

-0.000156776 

0.8 

22.8791 

0.0129460 

0.0329580 

-0.0009337 

-0.215893 

-0.000300122 

-0.000221710 

-0.000146949 

0.9 

22.8818 

0.0143893 

0.0311831 

-0.0063234 

-0.208231 

■0.000316451 

-9.000246179 

-0.000134615 

l.O 

22.8812 

0.0168825 

0.0289977 

-0.01.5,3842 

-0.197.574 

-0.00034.3217 

-0.000273061 

-0.000120122 

1.1 

22.8803 

0.020491’ 

0.0267376 

-0.0284821 

-0-18.5921 

-0.000.381016 

-0.000302055 

-0.000105483 

1.2 

22.8846 

0.0243030 

0.02.55578 

-0.0417119 

-0.179964 

-0.00042102.3 

-0.000.333949 

-9.75631e-05 

1.3 

22.8946 

0.025981- 

0.0267000 

-0.0449953 

-0.186824 

-0.000440711 

-0.000370519 

-0.000103697 

1.4 

22.8950 

0.0265152 

0.0275800 

-0.0428170 

-0.191549 

-0.000449334 

-0.000409584 

-0.000108446 

1.0 

22.8866 

0.027.5456 

0.0272887 

-0.0426774 

-0.18901 1 

-0.000462731 

-0.000447757 

-0.0001066.56 

1.6 

22.8752 

0.029396’ 

0.0263037 

-0.0461.528 

-0.182140 

-0.000484748 

-0.00048.3688 

-0.000101242 

1.7 

22.8640 

0.0319034 

0.0251308 

-0.0.524783 

-0.173920 

-0.000514.3.32 

-0.000517346 

-9.51958e-05 

1.8 

22.8542 

0.0347733 

0.0240901 

-0.060.3014 

-0.166228 

-0.000549028 

-0.000549196 

-9.0361 le-05 

1.9 

22.8459 

0.0376949 

0.0233346 

-0.0682171 

-0.159987 

-0.000586063 

-0.000579624 

-8.75.512e-05 

2.0 

22.8.381 

0.0403666 

0.0228641 

-0.0749269 

-0.15.5242 

-0.000622692 

-0.000608700 

-8.66960e-05 

00*^ 

22.7279 

0.0163541 

0.0.303799 

-0.0208336 

-0.151.376 

-0.000354967 

-0.0006.53371 

-0.000197299 


^ Cw = Co + CiWi + Co Wo + C3y^^ -h C4^^ + C.5wf + CeWjWo + C7W5. 

^ The oc represents the ternary system of NH3/HNO3/H0O (i.e., H2SO4 = 0). 


combinations (w^. W2; where 0.07c < wj + W2 < 85.0%) 
at cold temperatures, ranging from 180 to 270^K. The 
vapor pressures were then fitted into the following func- 
tions; 


0.0% < wi 4- W2 < 85.0%, ISO^K < T < 270‘M<, (3) 

B„(r. wn W 2 ) 


IiiPhno. == An(r, wi, w: 


InPn.o = Cw(r, wi, woj - 


T 

Dvv(r. wi, W 2 ) 


4) 


. o i 


where PuNOa PH2O are in mbar. Coefficients A^, 
B„, Cw and Dw were then parameterized into a mathe- 
matical function of wj and wj similar to that described 
by Luo et al. [1995]. The fitting results are simirnarized 
in Tables 1-4 for different r ratios, ranging from 0.0 and 
2.0 (r values are separated by an increment of 0.1). 

In Plate 1 the calculated vapor pn^ssures using e(|ua- 
tions (4) and (5) are compared with tlu>s(’ obtained from 
AIM2 for three r ratios (shown in Plates la, lb, and Ic). 
In general, PhsO calculated by r^quation (5) agrees well 
with the model of Clegg et al. [1998a|. However, the 


agreement of PhnOg using equation (4) is rather poor 
for most compositions. As shown in Plate 1. the cal- 
culated PHNO3 becomes increasingly inaccurate as the 
solution neutrality (i.e., r) increases. For some cases the 
difference between estimated PhnOs values and those 
predicted by AIM2 reached as high as a factor of 5. To 
ensure differences are minimized, Diff terms account- 
ing for the differences (deviations) of the “calculated" 
values (estimated from equations (4) and (5)) and the 
“rear’ values (obtained from AIM2) u'ere added to equa- 
tions (4) and (5): 


lnP[{NC)3 


A„(r, Wi, W2) -f 


Bn(r, w 

T 


Wj) 


A 


Diff„(r, vv'i , W2. T) (6) 


lnP[i,(; 


Cw(r. wi, W 2 ) 4- 


Dw(r. Wj, W 2 ) 

T 


+ 


Diffw (r, wi , wj, T). 


For a fix(^d i\ wi , and W2 the Diff correction terms ap- 
proximately follow simple polynomial functions in tem- 
perature: 

Diff„(r, Wi, wj. T) = rtTno + <7„iT + cXnjT’ (8) 
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Table 4. Dy, Coefficients for the Vapor Pressure of H >0^ 


r 

do 

di 

do 

d3 

d4 

d5 

de 

d- 

0.0 

-5850.24 

21.9744 

1.48745 

-44.5210 

59.6400 

-0.384362 

-0.6M671 

-0.208556 

0.1 

-5849.94 

19N594 

1.73322 

-39.7129 

58.0581 

-0.336545 

-0.606959 

-0.210361 

0.2 

-5849.15 

17,5776 

1.92716 

-34.3704 

56.5456 

-0.289409 

-0.568784 

-0.211815 

0.3 

-5848.04 

15.1511 

2.06677 

-28.5548 

55.1548 

-0.243024 

-0.530295 

-0.212915 

0.4 

-5846.79 

12.6262 

2.12829 

-22.4377 

54.0535 

-0.197733 

-0.491539 

-0.213544 

0.5 

-5845.48 

10.1120 

2.06175 

-16.4657 

53.5369 

-0.154503 

-0.452501 

-0.213397 

0.6 

-5843.72 

7 79449 

1.85460 

-11.3889 

53.6970 

-0.115030 

-0.41.).300 

-0.212396 

0.7 

-5840.63 

5.79437 

1.60810 

-7.6468 

53.9687 

-0.0803938 

-0.374279 

-0.211123 

0.8 

-5835.71 

4 .058 10 

1.46205 

-4.9098 

53.5707 

-0.0500282 

-0.335935 

-0.210378 

0.9 

-5829.33 

2.44674 

1.47151 

-2.5108 

52.2170 

-0.0225318 

-0.298645 

-0.210471 

1.0 

-5822.28 

0.S4270 

1.60302 

0.0831 

50.1179 

0.0033213 

-0.262577 

-0.211194 

l.l 

-5815.44 

-0."9036 

1.74377 

2.9989 

47.9281 

0.0280087 

-0.227587 

-0.211883 

1.2 

-5810.31 

-2.20859 

1.55496 

5.1396 

47.5949 

0.0492844 

-0.19.3207 

-0.210535 

1.3 

-5807.19 

-2.71615 

0.62592 

3.4557 

51.4599 

0.0605252 

-0.158668 

-0.204729 

1.4 

-5801.66 

-2.43122 

-0.39261 

- 1.5150 

55.7388 

0.0628630 

-0.124,306 

-0.198357 

1.5 

-5792.68 

- 1 . ?5557 

-1.06065 

-7.4818 

57.9126 

0.0612095 

-0.091134 

-0-193964 

1.6 

-5781.23 

-1.28487 

-1.31402 

-13.0971 

57.6578 

0.0586697 

-0.059876 

-0.191880 

1.7 

-5768.33 

-0.-8563 

-1.18923 

-17.6257 

55.2142 

0.0571507 

-0.0.30919 

-0.191856 

1.8 

-5754.83 

-0.”5269 

-0.74364 

-20.6765 

50.9331 

0.0578653 

-0.004421 

-0.19.3536 

1.9 

-5741.30 

-0.93290 

-0.03072 

-22.0791 

45.1436 

0.0615518 

0.019549 

-0.196612 

2.0 

-5728.09 

-1,44363 

0.91523 

-21.7913 

38.0640 

0.0686466 

0.040959 

-0-200889 


-5736.70 

-0.27462 

2.02241 

0.5073 

20.6699 

0.0337051 

-0.045872 

-0.179320 

^ D 

l-y W 

— do -h d; 

[ W 1 -f d‘2 W'2 -h d3 >%Wi + d4 ^ 

-h do 'W 

f + deWi w*_> 

+ dyW.]. 



^ The oc represents the ternary system of NH3/HNO3/H2O (i.e.. H0SO4 = 0 ). 


Ditfw(r. wi, T) = cTwo ^ ^wiT 4 - ciw-iT'^, (9) 

where a are second order polyiioiiiial coefficients fitted 
for a particular combination of w^, wj. and r. Tfie 
effects of including Diff terms in improving the vapor 
pressure fits are shown in Plates Id- If where HNO3 and 
H2O vapor pressures calculated from equations (6) and 
(7) (with polynomial coefficients of ecpiations (8) cind 
(9)) are compared against AIM2 results for the same r 
ratios examined above. As shown in Plates Id- If, Diff 
terms force the cahailated vapor pressuri's of both HjO 
and HNO3 to agree with AIM2 results. For r vahu‘S 
other than 0.0, 1.0. and 2.0, including Diff terms in 
\ai)or pressure relations product's a nearly exact agree- 
ment with AIM2 results. Thus, for A PCM we tabulated 
polynomial coefficients of a (lookup tables available as 
eh'ctronic supporting matt’riab ) for all possible combi- 
nations of wi. w>, and r. The w(‘ight [)er(*ents of W] and 


'Supporting lookup tables are available \ ia Web browser or via .-Vnonymous 
tTIM'rom Up: ko<mos agu oig. direelor\'“apend ’ (I semanie • ‘'anonymous , 
Password - 'guest'’); subdireetones in the Up site are arranged by paper 
number Inlonualion on searcliing and submitting eleetronic supplements is 
tbund at Imp www.agu.org pubs esupp about hlml. 


w *2 in the lookup tables cover and S])an the composition 
spectrum from 0.0% to 85.0% with 1% increments in 
composition. 

As temperatures cool, solubilities of trace gases, sucli 
as nitric acid, increase significantly, thereby dei^Ieting 
gas phase concentrations. Tlius it is essential to cal- 
culate the distribution of HNO3 between the gas and 
aerosol phases. To simulate the gas acuosol partition- 
ing of HNO3, ecpiations (6) and (7) are c'oiqdc'd with 
mass cons(‘rvation of HNO3: 

total HNO3 = Pji.NC);, + HN()3(j^^j (in rnol/in’^). (10) 

Three unknowns, (NH.| )rH2„rS()4^a(i) (or wi), HN()3(i^(,, 
(or wj), and PhnO:, - oan be uniciuely determined by 
equations (6). (7). and (10). For this work, the sarru' 
numerical scheme as that utilized in EQLTSOLX' II [Ja- 
cobson id ai. 199G] is applied to solve eiiuations (G)-(IO) 
iteratively. 

At iteration st(q>s wlu're W} and \\’> are not tabu- 
lated in th(‘ lookup tabh's, tlu' biliiu'ar (area w(‘ight(‘d) 
averaging [Jacobson, 1999a] is adopter 1 for iiit(U'pola- 
tion. For e.xainph', to estimat(* [)roperties at [wj. wjj = 
[25.3. 10.8], results ol)tained from four adjacent [)oiiits. 




Plate 1, Conij)arison of vapor pr('ssiir(‘s obtained from APCM parameterization and those 
predieted by A I M2 for tlire(‘ r ratios, (a, b, e) Vapor pressures ealculated by equations (4) and 
(5), and (d, f) r('sults (’alculat.ed i)y (‘cpiations (G) and (7). 
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Table 5. Composition Functions^ 






A 

B 

c 

D 






(SH^hSO^ 


0.00 < 

aw 

< 

0..30 






yi 



5,r22.^5602.33e-01 

9-6 193408.3360- 01 

l. 4720.57870 le+02 

-1.72291123670-02 


Vo 



4.5462384945er01 

-9. 187202S9I7O-01 

.3.3412155I66e+02 

■2.11246148250+02 

0.30 < 

aw 

< 

0.60 






>T 



9.7:j4-016855e-01 

-2.744425127.5e-f00 

-1.89665286140+01 

1.92120566470+01 


y-2 



1.014^647694e-f00 

-2.:t9l6961554e4-00 

-2.29879509860+01 

2.19200188650+01 

0.60 < 

aw 

< 

l.OO 






yi 



-5.375o389490e+02 

i. 35364.564.5.30-02 

-1.89512482810+01 

5.569210.5238e+02 


yo 



-4.29'2:431877e+02 

1.287.5856.3400-02 

-2.12759.335200+01 

4.50754787620+02 






NHi 



0.00 < 

aw 

< 

0.40 






yi 



1.3165236133e^02 

-1-01034156000+00 

-3.17268299680+01 

-1.40420757610+02 


Vo 



8.114y051028e+0l 

9.621338524.30-01 

1.05367666940+02 

-1.53609192270+02 

0.40 < 

aw 

< 

0.85 






yi 



8.57s2004669e-f0l 

-1.27956.359790+00 

2-99826614420+01 

-1.07762908850+02 


\ 2 



3.2723371194e^0l 

-1.32992178510+00 

-8.12459730750+00 

-2.45202617990+01 

0.85 < 

aw 

< 

1,00 






yi 



-1.4965599516e-01 

-2.02713825400^01 

-2.1065828892O+02 

2.06790766610+02 


Vo 



1. 663^062 187e^03 

-1.34654286080-01 

1.85922874440+02 

-1.84966766350+03 






;VH4//504 


0.00 < 

Uw 

< 

0.25 






yi 



8.9 14 !'5.56445o^0 1 

-9.22254497650 01 

2.16681577120+02 

-2.75523385310+02 


>'2 



5.984uS74514e-01 

9.37299480480 01 

2.72262245220+02 

-2.17215846220+02 

0.25 < 

aw 

< 

0,80 






y [ 



2. 8 75 34 59 72 5e-i- 00 

2-.55O7i561O90+OO 

-3.331686093.30-01 

-1.84942140120-03 


Vo 



3.555o024880e-00 

-2.0910350.5000+00 

-1.0492539829O+01 

7.89230135130+00 

0.80 < 

aw 

< 

1.00 






yi 



-1.422'3.39985o-00 

2.927723 14320+01 

-1.721.50580880+01 

1.84416626360+01 


>'2 



-7.04-3234774e-0l 

2.08.32745270O+01 

-2. .32895939290+01 

2.39011375590+01 


'' y - AaJ,^ -h CRw + D, Read 5. 1229560‘233e+01 as 5.1229560233 x 10^ All parameterizatioiis are 
valid for 190" K < T < 2G0" K only. Molality rn is calculated by iii(aw,T) — yi(aw) + (T - 190) 
” //i(fLr)] /~0> where a^ is water activity (relative humidity expressed in fraction). Compo- 
sition functions for H>S 04 and HNO 3 in the same mathematical form are given by Tabazadeh et at. 
[1997a. 1997b], 


[25.0. 10.0]. [25.0. 11.0), [26.0, lU.U]. [26.0. 11.0], are 
weighted and a\'eraged, de[)eiuliiig (;ii liow far the de- 
sired point is located away from the four hxed points. 
Similarly, a simple linear interpolation (t>f avc^raging 
prcflitnions from two adjacxmt r ratios) is implemented 
if initial r falls inhetween the 0.1 increnuMUs (.see Ta- 
ble IJ. Siric(‘ the intervals are small (P/{ for w’l and wj 
and 0.1 for r), the two interpolati(m methods yield rea- 
sonable ri'sults that are in close agreement witli AIM2 
in predicting HNO3 and HjO vapor pressures ov('r the 
N H;j / H 2 S C) 1 / H N 0 3 / H-2 O system . 

For conditions where HjSO.j mixing ratio is Z(‘ro w(' 
defin(‘ w^ and W2 as the weight ptTcents of NH,jN()3{.^,jj 
and HN03|^m). resptx'tivtdy, and rf jn^at the above st(*ps 
to obtain vapor pressures and s(diiri(jii coin])ositions for 
the NH3/HN03/H-20 systcMii. 


3.2. An Update of EQUISOLV II 

In addition to APCM parameterization an equilib- 
rium model basixl on the binary activity approach [Stel- 
son and Seinfeld. 1981: Pilmi.s aiid Seinfeld. 1987: ll exl- 
er and Seuifeld. 1991; ATm and Seinfeld, 1995: Jacob- 
son et ai. 1996; Jacobson. 1999b], EQUISOL\’ II [Ja- 
cobson et ai. 1996; Jacobson. 1999b]. is updated for 
application at colder temperatures. Tlie updated EQ- 
I'lSOLV II will be compared against AIM2 and APCM. 
The important (H[uilibrium ecpiations to solve in EQLT- 
SOLV II for the NH3/H>S().,/HN()3/H20 system are 


HSO 




Hr +so:r r 

(a(|; ‘Ufup 


:ii) 


H-N-0:hk, + NO^; 


3(a(p ' 


( 12 ) 
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Tabic 6. Moan Binary Activity Coefficients^ 


QkiT) = 3o + 0iT + 

Coefficients 

/^o 

iii 




H { n = 5 ) 



1 . 1 55060753 le-|-01 

-8.6598065340O-02 

1.61003351326-04 


2.60456^6534e-f01 

1.85701 18132c -01 

3.40.329316416-04 

a 2 

9.6192566768e-i-00 

-6.43141207286-02 

1.1363959743e-04 


-1.4330104905e-r00 

9.34704903936-03 

-1.632107.39286-05 

Odj 

9.5181767818e-02 

-6.1 140131 1740-04 

1. 06058569 19e-06 


-2.3568483179e-03 

1.493602 1604O-05 

2.5755879088O-08 



( ) 2 ‘8334 ( ?i = 3 ) 



3.4428:516946e-|-00 

-4..53283216546 02 

1.02021375296-04 


-1.00306.30804e-t-00 

6.0339372495C-03 

-1.2162329645c 05 

a 2 

4.2925175339e-02 

-3.0262230896O-04 

5.9904 1767680-07 


-6.88 1.3265363c 04 

5.003 16887980-06 

'9. 746025858 7e-09 



O 3 { 71 = 4 ) 


ao 

2.7037375056e4-00 

- 1.78H0171498O-02 

2.84977783930-05 


4.y424685844e+00 

2.16218493340-02 

- 2.2789057453e-05 


9.0584578 176e-0i 

-4.29471613830-03 

4.79719055I6O-06 


7.755469I781e-02 

3,76 1 8465057e -04 

4.3345626096O-07 

Qi 

2.4783745998e-03 

1 .2123 1 65585o-05 

1.41814244910 08 



\ H.x H S(J.^ {n = ^j 



9.960 14774510 01 

-1.12969512910 -02 

1.74527 1 67440 -05 

0^1 

1.484583 145 le+00 

-9.3756384679O-03 

I.8989492030O 05 

0-2 

2.9273249432e-01 

1.3469653 1740-03 

-2.27189173596-06 


2.19885859490-02 

-8.1477260064O-05 

1 .04703568716-07 

a.i 

-7.53971124320-04 

2.209686396.3O-06 

-1.60518210786-09 


9.7902066682e~06 

-2.15603962210 08 

-3.3248.336670e-12 


ak(T)rn^/“, Read 1.1550607531e-h01 as 1.1550607531 x10^. All 


Ai. vyiiiy. luu rtiiiviLv lui me 

HjS 0.,/H20 s\stem \va.s not parameterized but instead tabulated in the computer 
coilos. 


NH, 




+ H+ , 
(an) 



(13) 


where Kj is the etjiiililiriuin constant and (g) and (aq) 
refer to gas- and atpieous-phase sjiecies. respectively. 
The first dissociation step of Hj.SO.i (H2SO,,||, — >• 
assumed to bf- complete. 

Temp('raturc-dei)(>ndent water and solute activity co- 
efficients in EQUISOLX' II are modified at cold tem- 
peratures (190'’K-26U"K) by parameterizing data from 
the model of Cleyg et al. [1998aj. Tin* ttmipmature- 
dependent water activities of five elc-ctrolytc's involved 
in etiuations (11)-(13), (H-jSO.,. HX():j. (.\H.,)2SO.|. 
XH,.\();j, XH.jHSO.,) are improved in EQITSOLV II. 
For consistency with iirevious work by Tabazadc.h c.t al. 
[199ia. 1997b) the binary solution ( (mqiositions are pa- 
ranieterizi'd into 


= yi(avv) -b ^ 

70 


1 ) 


where in is the molality of the binary electrolytes, a^ is 
the water activity (relative humidity e.xpressed as frac- 
tion), and y is a function of a« only (see Table 5). 
Mean binary activity coefficients (yiy) were parame- 
terized by converting mole fraction actic ity coefficients 
(f) of Clegg et al.'s [1998a] model into molality base 
(7) values according to eciuation (15). and using = 

[Pitzer. 1991); 

(o) 

wlicrc and arc tlu* stoichionuaric co(dficuuits of 
the binary electrolytes (for example, = 2 and i/_ = 1 
for (NH.| )jS04), Mw i-s tlu‘ nujlecular weight of wat(u\ 
and sumnitition is over all solute species. The nu'an 
activity (‘oefficients vv(‘re then fitted into sirriph' poly- 
nomial functions (given in Table 6) of temperatuia* and 
molality (where and .4 are polynomial coc'fficients): 
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Tabl e 7 . Model Conditions^ 

NH3 


Cases 


Background 

Volcanic 

Background 


Polluted 


Background 


H2SO4 


r = 0.0 
r = 0.0 

r = 0.0 
r = 1.0 
r = 2.0 
r = 0.0 
r = 1.0 
r = 2.0 

r 0.0 - 2.0 
r = OC 


Laaksonen et al., 1997]. 


k = 0 

ak(T) = /?o + ^iT -r • 


NH3, ppt 

H2S04. ppt 

HN03, ppt 

H2O, ppm 

Illustration 

Stratospheric Simulation (Ptotai -^0 mbar ) 



0.0 

500.0 

10,000.0 

5.0 


Figure la and lb 

0.0 

20,000.0 

10,000.0 

5.0 


Figure 2a and 2b 

Upper Tropospheric Simulation {Ptotai mbar ) 



0.0 

100.0 

100.0 

5.0, 50.0, 

500.0 

Figure 3a 

100.0 

100.0 

100.0 

5.0, 50.0, 

500.0 

Figure 3b 

200.0 

100.0 

100.0 

5.0, 50.0, 

500.0 

Figure 3c 

0.0 

100.0 

2,000.0 

5.0, 50.0, 

500.0 

Figure 4a 

100.0 

100.0 

2,000.0 

5.0, 50.0, 

500.0 

Figure 4b 

200.0 

100.0 

2,000.0 

5.0, 50.0, 

500.0 

Figure 4c 

HNO, Uptake (Ptotai =^00 mbar ) 




0.0 - 200.0 

100.0 

100.0 

5.0, 50.0, 

.500-0 

Figure 5 

100.0 

0.0 

200.0 

5.0, 50.0, 

,500.0 


pressure, and the mixing ratios used are atmospheric observations [Tabazadeh et ai, 1998; 

r)m‘'''^ 

( 16 a) 

The above parameterizations 

induce relative errors 

of no more than a few percents 

for both water and so- 


lUte aCLlVlLies. J.11C 

( 16 b) ity regions where ion activities usually exhibit a com- 



Fieure 1. Comparison of various model predictions of (a) vveiglit percents and (b) aerosol 
volume under background conditions in the stratosphere (see Table 7 for model conditions). 
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7J C - 


2:1.98 J?.28 26.9^ 19 t2 
Cci’''"espo'”2 ' 2 


'-i.39 10.52 /.9C 5.31 

= H (%) 


/J.c/ 51.98 J.^28 26.9*i '9 62 ' 4^9 '0 62 90 5.9* 

Correspond n.g RH { 7 ,) 


Figure 2. Comparison of various model predictions of (a) weight percents and (b) aerosol 
volume under volcanic conditions in the stratosphere (see Table 7 for model conditions) . 


plicated behavior. In addition to the improvement of 
binary water and solute activity coefficients, equilib- 
rium constants of equations (ll)-!l3) were obtained 
from Clegg et al [1998a] and converted into suit- 
able rnolality-based units for use in EQLTSOL\' II. 
These modifications ensure perfect agreement of EQUT 
SOLV II and AIM 2 at a binary level where electrolyte 
solutions are made of H>0 plus only one electrolyte. 

4. Model Intercomparison and 
Evaluation 

4.1. H.JSO 4 /HNO 3 /H 2 O 

For the stratospheric system of H 2 SO 4 /HNO 3 /H 2 O, 
the compositions obtained from APCM are compared 
with previous formulations [Carslaw et al., 1995b; Lilo 
et at., 1995: Tahazadeh et ai. 1994j under background 
and volcanic stares (see Table 7 for model conditions 
used). Results from model intercornparisons between 
six different ternary aerosol models are shown in Fig- 
ures 1 and 2. For the background stratosphere, APCM 
with Diff correction terms included is in good agree- 
ment with AIM2 ^Car'.Hlaw et ai. 1995a]. The predicted 
weight percents by APCM overlap with AIM2 almost 
exactly throughout the entire temperature range shown 
in Figure la. whert?as previous models deviate from 


ALM2 either at the initial, in the middle, or at the 
final stage of aerosol growth (Figure lb). Under vol- 
canic conditions (H 0 SO 4 increases from 0.5 to 20ppb), 
all models produce nearly identical results with only 
slight variations (Figures 2a and 2b). 

4.2. XH3/H2S04/HN03/HoO 

In the upper troposphere, aerosol compositions are 
examined for both background and polluted states (see 
Table 7 for model conditions). Under polluted con- 
ditions, HX 0;3 mixing ratios are elevated well beyond 
background values of lOOppt [Laaksonen et ai, 1997], 
mainly because of com'ective transport of polluted boun- 
dary layer air dirt‘ctly into the upper troposphere. \ ari- 
ations in the aerosol composition predicted by three 
models, APCM. modified EQlTSOL\’ II (with the im- 
proved activity data), and AIM2 [Clegg et al., 1998a], 
are illustrated in Figures 3 and 4. The points of ice sat- 
uration are marked ms vertical dotted lines in Figures 3 
and 4. For all cases studied and compared, APCM 
yields compositions that are in close agreement with 
AIM2. The modified EQITSOLV II with new activ- 
ity data portonns better in the regions where relative 
humidity (RH) is high. The deviations at lower rela- 
tive humidit}' are lik(4y caused by the unct^rtainties in 
mixed activity coefficients calculated from simple mix- 
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Figure 3. Variation of aerosol compositions as a function of temperature (or relative humidity) 
under background conditions in the upper troposphere (see Table 7 for model conditions), (a) 
r = 0 . 0 , (b) r = 1 . 0 , and (c) r = 2 . 0 . 


irig rules [Dioinley, 1973: Kusik and M€is.iuc7\ 19/8] in 
EQUISOLX' II. 

We hav(‘ also compared predicted niolalitit^s of indi- 
vidual com[)onents, NH^, H 2 SO 4 (= HSC)^ -h S();j ). 
and NO 3 , in solution (instead of weight percents of 
vvi and wj) and found that the differences between 
APCM and AIM2 are no more than a few p(?rcent ex- 
cept at the very low or high solute concentration re- 


gions. In other w'ords, predictions by APCM are less 
accurate only in regions where weight percents approacli 
the lower or upper limits (P^f and 85%). Nevertheless, 
for very low w'eight percents, concentrations are too di- 
lute to be significant, whereas for very high weight per- 
cents, relative humidities are sufficiently low that am- 
moniated salts, such as (NH 4 ) 2 S 04 or (NH 4 ) 3 H(S 04)2 
(letovicite), would most likely precipitate in solution 
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(%) for ^20 = 52pm (%) for H20 = 50ppm 

99.5 67.1 45,3 3\6 22.0 96.3 /8.9 57.1 41.6 30,6 

I 1 I ^ ^ 1 1 


RH (%) for H20 = 500ppm 
96.6 73.7 56.6 45.7 
( 1 ^ 1 



Figure 4. Variation of aerosol coinposit.ions as a function of temperature (or relative humidity) 
und(T polluted conditions in the upper troposphere (s<’e Table 7 for model conditions), (a) 
r = 0.0, (b) r = 1.0, and (c) r = 2.0. 


[Tahazadeh and Toon, 1998]. For (example, crystalliza- 
tion of (NH 4 ) 2 S 04 occurs at -- 35% RH at room tem- 
perature [A a et ai, 1998). At colder temperatures, crys- 
tallization of (XH 4 ) 2 S ()4 will probably o(‘(“ur at slightly 
higlun* RHs. In addition, for v(tv concentrated solu- 
tions (W[ -h Wo > 85%), aerosol compositions predicted 
by AIM2 are often outside the rang(‘ of model validatioii 


[Clf'(j(j ct al., 1998a]. For the reasons mentioned above, 
the limits imposed on the weight percents in the APCM 
are not a seritnis drawback. Also, as sliown in Figures 3 
and 4, the APCM compositions are most accurate near 
regions where ice reaches saturation in the atmosphere. 

In addition to aerosol composition the APCM can be 
used to calculate the extent of HXOj uptake by up- 
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Figure 5. Uptake of H NO 3 in sulfate-based aerosols as a function of sulfate neutralization 
under background conditions in the upper troposphere (see Table 7 for model conditions), (a) 
H.O = 5ppnn (b) H2O = oOppm. and (c) HjO = 500ppm. (Dashed lines represent the uptake 
by pure XH.1XO3 aerosols (i.e., H2SO4 ~ ())) 


per tropospheric aerosols. Aircraft field experiments 
hav'e shown that ammoniated particles are abundant 
in the upper troposphere [Talbot et a/., 1996, 1998; 
Tabazadth et ai. 1998]. Here vve examine how HNO3 
u])takc may be affected by the presence of ammoni- 
ated aerosols in the upper troposphere. For simulations 
the H2SO.1 and HXO3 mixing ratios are fixed at back- 
ground levels of 100 ppt each. XH3 concentrations in 
the APCM are varied by increasing (or decreasing) r 
to account for changes in solution neutrality (or pH), 
ranging from pure H0SO4 solution droplets (r = 0.0) to 
fully ammoniated systems consisting of only (NH4)2S04 
(r = 2.0) (see Table 7 for model conditions used). 


In Figure 5. the uptake of HXO3, expressed as a frac- 
tion (i.e.. ratio of concentration in the liquid phase to 
the total initial concentration), is shown for three differ- 
ent assumed w-ater vapor pressure profiles. As expected, 
partitioning of HNO3 in sulfate-based aerosols depends 
strongly on solution neutrality. In general, greater up- 
take occurs at lower temperatures and higher relative 
humidities. For example, at temperatures lower than 
210^K (which corresponds to a relative humidity of 78% 
for Ptotai = 200mb and H2O = 50ppin), a significant 
fraction (> 90%) of HNO3 resides in fully neutralized 
ammoniated solutions, compared to only < 20% in the 
pure sulfate system (Figure 5b). It has been shown that 
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high levels of HNO3 in solution may cause precipita- 
tion of amnioniated and/or nitrated salts, which could 
change the mode of ice formation from homogeneous to 
heterogeneous nucleation [Tabazadeh and Toon, 1998 ]. 

5. Conclusions 

The ion interaction model of Cltgg et ai [ 1998 a] for 
the system of NH3/H2SO4/HNO3/H2O has been pa- 
rameterized into a compact model well suited for in- 
corporation into large-scale atmospheric models. The 
aerosol physical chemistry model (APCM) reproduces 
the AIM2 results of Clegg et ai • 1998 a] for a wdde 
range of conditions in the upper trr>posphere and lower 
stratosphere. Model intercomparisons show that for the 
ternary system of H2SO4/HNO3/H2O solution compo- 
sitions obtained from APCM are in better agreement 
with those obtained from AIM 2 than previous formula- 
tions. For the quaternary system of NH3/H0SO4/HNO3 
/H2O, APCM results are also in good agreement with 
AIM2 predictions, particularly near the regions of ice 
saturation where the influence of amrnoniated particles 
on the ice nucleation process is of interest. Extension 
of APCM to include other features such as calculations 
of deliquescence relative humidity and precipitation of 
solids in solution are under way. 
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Abstract. 

Satellite observations of water vapor and aerosol extinction along with temperature 
trajectory calculations are analyzed for the Southern Hemisphere winter of 1992 
in order to determine the onset, extent and duration of dehydration within the 
polar vortex. Our investigation utilizes measurements of water vapor from the 
Microwave Limb Sounder (MLS) and extinction from the Cryogenic Limb Array 
Etalon Spectrometer (CLAES) both onboard the Upper Atmosphere Research 
Satellite (UARS). Evidence for dehydration is seen on potential temperature 
surfaces from 420-520 K (approximately 16 to 22 km). The horizontal extent of 
the dehydrated area at 465 K encompasses up to 35% of the total vortex area. 
Based on MLS observations and temperature statistics, the onset of dehydration 
occurs between late June and early July, and dehydrated regions persist well into 
November. A comparison of CLAES aerosol extinction measurements and model 
calculations of extinction suggests an average ice particle number concentration 
and size of 10~“-10~^ cm“^ and 10-30 /im, respectively. We show that the 
difference in timing of the onset of dehydration found here an in a recent analysis 
of Polar Ozone and Aerosol Measurement III (POAM) observations is due to the 
latitudinal sampling pattern of the POAM instrument. 


1. Introduction 

Polar stratospheric clouds (PSCs) form in the low temper- 
atures of the winter polar vortex and have a direct influence 
t)n atmospheric chemistry by their involvement in heteroge- 
neous ozone chemistry [Solomon, 1999). PSCs also act as 
a sink for gas phase water vapor and nitric acid resulting in 
stratospheric dehydration (e.g., Kelly et ciL, 1989; Hofmann 


and Deshler, 1991; Nedoluha et ai, 2000] and denitrifica- 
tion [e.g.. Toon et ai, 1986; Santee et ai, 1999; Tabazadeh 
et 2000a]. 

Understanding the physical properties of the clouds and 
the meteorological conditions under which they form is im- 
portant for studies of the polar vortex. Particle size and com- 
position influence the optical properties of the clouds, which 
affect remote sensing [Toon et al., 1990; Browell et ai. 
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1990]; composition of the aerosols impacts constituent bud- 
gets and chemical reactions [Solomon. 1999]; and chemical 
processes which are dependent upon relative humidity will 
be affected by the deficit of water vapor due to dehydration 
[Hanson et ai, 1994]. PSCs are separated into two types 
[WMO, 1999]. Type II PSCs are composed of water ice and 
form at temperatures below the ice frost point. Type I clouds 
form at temperatures above the ice frost point and consist of 
both liquid and solid nitric acid-containing cloud particles. 
We focus our study on the atmospheric dehydration process 
which occurs as a result of ice particle growth and sedimen- 
tation, thereby removing water from the atmosphere. 

The Airborne Antarctic Ozone Experiment (AAOE) in 
August-September 1987 showed the first evidence of strato- 
spheric dehydration in the Antarctic vortex [Kelly et al. 
1989]. Balloon measurements in the Antarctic have also 
shown profiles with water vapor depletion [Hofinann et ai. 
\99i \ Rosen et ai, 1988; Vdmel etai. 1995]. Using balloon 
measurements over McMurdo in 1994. Vdmel et ai [1995] 
conclude that dehydration starts around mid June and per- 
sists into November. Overall, only a limited number of in 
situ measurements have been made o\er the Antarctic re- 
gion over the course of several winters. Satellite observa- 
tions can provide a more extensive view of the polar re- 
gions. Springtime dehydration has been noted in the analysis 
of Halogen Occultation Experiment (H.ALOE) data [Pierce 
et a/., 1994; Rosenlof et ai, 1997]. Santee et al. [1995] re- 
port on the interhemispheric differences in wintertime water 
vapor from the Microwave Limb Sounder (MLS) and Ricaud 
et ai [1995] deduce PSCs from MLS w ater vapor and Cryo- 
genic Limb Array Etalon Spectrometer (CLAES) extinc- 
tion measurements from several days in August-September 
1992. Evidence for wintertime dehydration is also found in 
measurements from the Atmospheric Trace Molecule Spec- 
troscopy (ATMOS) instrument [Manner; et ai, 1999], and 
Nedoluha et ai [2000] report on Polar Ozone and Aerosol 
Measurement III (POAM) observations which show dehy- 
dration occurring in the Southern Hemisphere (SH) winter 
of 1998. 

In our study, we examine satellite measurements of strato- 
spheric water vapor and extinction in the SH winter polar 
regions to identify the onset, extent and duration of dehydra- 
tion. We correlate water vapor concentrations w ith aerosol 
extinction and use temperature trajectories to examine where 
and w hen the dehydration process occurs. We compare mea- 
sured and calculated extinction to infer ice particle size and 
number density in the dehydration column. 


2. Observations 

Water vapor data for this study come from the MLS in- 
strument onboard the Upper Atmosphere Research Satellite 
(UARS). We use the prototype water vapor product (version 
104) described by Pumphrey [1999], which gives water va- 
por mixing ratios in the stratosphere from 0.1 to 100 hPa 
with a vertical resolution of 3-4 km. Although this prod- 
uct is not an official UARS project MLS product, its validity 
and usefulness in scientific studies has been demonstrated 
[Pumphrey et ai. 2000; Pumphrey. 1999; Manney et ai. 
1998]. The water vapor measurements are available from 
September 1991 to April 1993, which includes only one SH 
winter. 

The CLAES instrument, also onboard UARS, measures 
aerosol extinction coefficients [Mergenthaler et ai. 1997] 
which have been used to identify PSCs [Mergenthaler et ai. 
1997; Ricaud et ai. 1995; Massie et ai. 1994] and to de- 
termine volume and surface area densities of type I PSCs 
[Massie et ai. 1997, 1998]. In identifying PSCs, extinction 
values are generally used qualitatively and threshold values 
for PSCs, depending on pressure, are determined relative to 
temperature and pre- winter observations. No information on 
particle size and concentration is given with the aerosol ex- 
tinction data product. CLAES was operational from Septem- 
ber 1991 until May 1993. Our analysis uses the CLAES ver- 
sion 8 extinction coefficients at 780 cm~ ^ 

MLS and CLAES measurements are nearly coincident 
in time and space, each with a horizontal resolution of 
~400 km. The vertical profiles of water vapor and extinc- 
tion have been interpolated from the standard UARS grid to 
potential temperature surfaces using United Kingdom Mete- 
orological Office (UKMO) temperatures. UKMO temper- 
atures are also used for the analysis in sections 2 and 3. 
Comparisons of the UKMO temperatures with radiosonde 
observations show that systematic biases are less than 1 K 
throughout the SH winter [Manney et ai. 1996]. 

We focus on the UARS SH view ing periods for one year, 
starting in April 1992 ( 1-30 April 1992, 2 June-12 July 1992, 
14 August-20 September 1992, 30 October-28 November 
1992, 10 January-8 February 1993, and 20 March- 1 April 
1993). We note that there are significant data gaps in the 
June-July viewing period due to problems with the UARS 
solar array which required the MLS and CLAES instruments 
to be turned off. Data availability limits this study to only 
one SH winter and we therefore do not assess year-to-year 
variability of dehydration. Although the SH winter of 1992 
was not necessarily a typical year since there were high lev- 
els of stratospheric aerosol from Pmatubo, the duration of 
cold temperatures and minimum temperatures were typical 
of wintertime conditions. 
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The water vapor annual cycles on potential temperature 
surfaces, compiled from the high-latitude MLS data, are 
shown in Figure 1. Daily averages are from data in the 
70°S-80®S latitude range for one year starting in April 1992. 
The 420 K potential temperature level is below 100 hPa in 
January through March, resulting in a lack of observations 
on that surface during that lime. Potential temperature sur- 
faces from 420 to 585 K show a minimum in water vapor 
in the winter months. These potential temperature surfaces 
cover pressure levels of 100 to 15 hPa. Since the water vapor 
data set does not e.xtend below 100 hPa. we cannot observe 
the lower altitude limit of the dehydration. Nedoluha et ai 
[2000] show evidence for dehydration down to 12 km. At 
the higher potential temperature surfaces, the maximum in 
the annual cycle is found in late winter (.August-September). 
The depleted winter values at and below 585 K are a result 
ot cold temperatures found within the polar vortex. Average 
temperatures in July on the 465 K potential temperature sur- 
face are near 185 K. below the ice frost point (188 K). Ice 
clouds form at these cold temperatures and deplete the gas 
phase water. 

The daily averages of aerosol extinction from CLAES are 
also shown in Figure 1. In general, they are anti-correlated 
wdth the w'ater vapor in the winter months from 420 to 585 K. 
High extinction values are seen from mid June through late 
August when the water vapor mixing ratios are low, rein- 
forcing the observation that the deficit in gas phase water 
coincides with the presence of PSCs. Extinction values in 
September are lower than the pre- winter values. This fea- 
ture in the annual development of extinction has been seen 
by other (.ibservations and is due to the cleansing of the at- 
mosphere by the subsidence of PSCs and by descent in the 
vortex [Kent et ai, 1985; Thomason and Poole, 1993; Ran- 
dall et ai, 1996; Nedoluha et a!., 2000] 

Due to limited MLS water vapor measurements in late 
June and early July, we cannot observe the precise onset of 
dehydration. The measurements indicate that dehydration 
begins sometime after June 18. Extinction values are already 
high prior to June 18 due to Type I PSCs [Tabazadeh et ai, 
2000a], which form at temperatures higher than the ice frost 
point. A deficit in the water vapor mixing ratio is clearly 
seen in the July 10-12 measurements for the 420 to 585 K 
surfaces. At 465 K. daily average values have dropped from 
4 ppmv in mid June to 2.7 ppmv on July 10, indicating that 
the onset of substantial dehydration occurred sometime af- 
ter mid June and before 10 July. Water \apor values at the 
420 K and 465 K levels drop further between mid July and 
early August when the next observation period begins. In 
August-September, daily mean values are between 2.3 and 
2.8 ppmv at 465 K. Extinction values are still elevated in mid 
August but have diminished by early September. Measure- 


ments on the 585 K surface show a drop in the water vapor 
and an elevation of extinction in mid July, but values have 
recovered by August, indicating that the initial drop in water 
vapor seen at this level was reversible or that mixing ratios 
are replenished by the descent of moister air from above. 

Figure 1 depicts only daily averages in the 70°S-80°S lat- 
itude range. Additional information can be found in the dis- 
tributions of water vapor and extinction measurements. Fig- 
ures 2a and 2b show distributions of water vapor and aerosol 
extinction in the region of 50''S-80®S for different temper- 
ature bins (T<185 K, 185<T<188 K, 188<T<195 K, and 
1 95 <T<2 10 K (points w ith temperatures greater than 210 K 
are excluded). We also show the distributions of water vapor 
and extinction for the 70'^S-80°S latitude bin (for all temper- 
atures). The T< 1 85 K bin represents air below the ice nucle- 
ation point (typically 3-4 K below the frost point [Tabazadeh 
et ai, 2000b; Chang et ai, 1999]) and the 185<T<188 K 
bin represents observations between the frost point and the 
nucleation point. The distributions for data taken between 
.April and November 1992 are separated according to the 
CARS yaw cycle. 

We use the high latitude April measurements, represent- 
ing pre-wintertime observations, to define low water vapor 
and high extinction values later in the paper. The water vapor 
distribution peaks at 4.3 ppmv in April and all the extinction 
values in this month are less than 0.0009 km“^ In June- 
July, the mode of the water vapor distribution for T<185 K 
is 2.3 ppmv and the 185<T<188 K observations have two 
peaks at 2.3 ppmv and near 4 ppmv. For T> 1 88 K the distri- 
butions are centered at 4 ppmv and fewer than \07c of the ob- 
servations are less than 3.5 ppmv, indicating that by this time 
the dehydrated area has not yet spread beyond the stability 
region of ice particles. Therefore, we define dehydrated ar- 
eas to have water vapor mixing ratios less than 3.0 ppmv 
based on the April measurements, and “high" aerosol ex- 
tinction is defined to be values larger than 0.001 km~^ The 
aerosol extinction threshold is slightly higher than the value 
used to identify PSCs at 46 hPa (0.00075 km“^) by Mergen- 
thaleret ai [1997]. Extinction values over 0.002 km~* have 
large uncertainties due to the high optical depth of the PSCs 
along the limb viewing path [Mergenthaler et ai, 1997]. 
However, the uncertainties do not affect the classification of 
‘high" aerosol amount. The June-July distributions show 
that nearly all the observations which coincide w ith temper- 
atures below 185 K are dehydrated. 25^ of the observa- 
tions with 185<T< 188 K are dehydrated. For temperatures 
above 1 88 K, there is no evidence of significant depletion of 
the water vapor. 

The extinction distributions for the June-July observation 
period show’ evidence of both Type I and Type II PSCs. 
The temperature bins with T<195 K all show observa- 
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tions of extinction above the April distributions. For the 
185<T<188 K temperature bin, 75% of the observations 
have extinctions larger than the pre-winter values. This is 
significantly larger than the 25% of water vapor measure- 
ments that show a deficit in this temperature bin, suggesting 
that type I PSCs are contributing to the elevated extinction. 
The coldest temperatures coincide with the highest extinc- 
tion values. All of the observations with T> 195 K show low 
extinction values, less than 0.001 km“*. 

By the August-September viewing period, the distribu- 
tions of water vapor mixing ratios for all temperature bins 
with T<195 K peak near 2.2 ppmv. For 195<T<210 K, the 
majority of the observations are centered around 4 ppmv but 
there are measurements with mixing ratios less than 3 ppmv 
in this temperature range as well. The extinction distribu- 
tions during this time have a small percentage of observa- 
tions with extinctions above 0.001 km"’ and these coincide 
with the coldest temperatures (T< 188 K). Thus by this lime 
in the winter, there exist dehydrated mixing ratios that do 
not directly coincide with low temperatures or high extinc- 
tion values. 

Although temperatures at high latitudes are warmer by 
November, water vapor mixing ratios have recovered only 
slightly from the dehydration. The distribution for TO'^S- 
80°S peaks at 3.0 ppmv. While some water vapor recovery 
has taken place in November, the overall mixing ratios are 
still substantially below the April mode value (~4.3 ppmv) 
before any dehydration occurred. 

3. Extent and Timing of Dehydration 

In this section we quantify the horizontal extent of the 
dehydration, its development throughout the winter months, 
and its relationship to the areas of cold temperatures, high 
extinction, and the vortex. To do this, we calculate the area 
encompassing the vortex and the regions of dehydrated air, 
high extinction, and cold temperatures. Criteria for dehydra- 
tion (mixing ratios less than 3.0 ppmv) and “high’’ extinction 
(larger than 0.(X)1 km"’) are determined from the frequency 
distributions in the early winter months and before the onset 
of winter (Figures 2a and 2b). The water vapor distributions 
for early winter show two peaks, separated near 3 ppmv, and 
April measurements all exceed 3 ppmv (Figure 2a). The ex- 
tinction distributions in early winter show a distinct separa- 
tion at 0.001 km"’ and pre-winter values are all below that 
threshold (Figure 2b). The cold temperature criteria are the 
frost point of ice (T<188 K) and the ice nucleation tem- 
perature (typically 3 K less than the frost point, T<185 K 
[7(ihazadeli ct ai, 2000b; Chang et al., 19991). We use a 
scaled potential vorticity value of 1.4 s"’ to approximate the 
vortex boundary. F^igure 3 shows the calculated daily areas 


from April through November on the 465 K potential tem- 
perature surface. 

The area of temperatures below the frost point of ice be- 
gins to develop in mid June. At this time high extinction 
values are present but there is no evidence of dehydration . 
The area of high extinction values drops on June 24 and then 
increases again. Tabazadeh et ai [2000a] show that deni- 
trification at 450 K occurred over 12 days from June 12-24. 
Temperatures below the ice nucleation point are present over 
a small area in mid June but do not persist or develop over 
larger areas until after June 27. By July 10, the area of dehy- 
dration coincides roughly with the area of temperatures be- 
low the ice nucleation point and is similar to the area of high 
extinction. When water vapor measurements are available 
again in mid August, the dehydrated area exceeds the cold 
temperature area and persists at 12 x 10^ km^ through most 
of September, consistent with observations by [Vdmel et ai, 
1995] in which they find that the water vapor profiles from 
early August and early October in 1994 at McMurdo Station 
are similar. From mid August through late September the de- 
hydrated area is roughly 30%-35% of the total vortex area. 
Cold temperature areas decrease rapidly at the end of Au- 
gust. The area of high extinction is much less than the dehy- 
drated area during August and September, indicating that the 
dehydration persists after the ice clouds have precipitated. 
The dehydrated areas persist and remain at nearly a constant 
value (from early August to late September) after the area of 
high extinction (indicating the presence of clouds) has van- 
ished in early September and after the vortex area begins to 
recede in early September. The dehydrated area decreases 
rapidly after November 15, and the water vapor values start 
a rapid recovery. The water vapor mixing ratios at 465 K 
(Figure 1) show an increase between mid September and the 
beginning of November similar to what is observed for nitric 
acid [Santee et ai, 1999]. The dehydrated area in Figure 3 
does not show this recovery occurring since the area will not 
start to decrease until mixing ratios surpass 3 ppmv, indicat- 
ing that mixing ratios are attaining values similar to those 
present before temperatures fell below the ice frost point. 

4. Temperature Statistics 

The interpretation of the time development of the area 
of dehydration compared to the area of cold temperatures 
(Figure 3) is limited in that it does not include the temper- 
ature history of air parcels. Tabazadeh et ai [2000a] uti- 
lize the concept of 'TSC lifetimes" to examine how long an 
air mass persists at temperatures less than the condensation 
points of nitric acid trihydrate (NAT), nitric acid dihydrate 
(NAD), and water ice. In this section, we use temperature 
statistics from National Centers for Environmental Predic- 
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tion (NCEP) trajectory calculations to look at the onset of 
temperatures which meet the ice nucleation conditions for 
the formation and growth of water ice clouds. 

For each day in June through September, 20 points were 
homogeneously distributed within the 188 K temperature 
contour at 50 hPa. Twenty forward and backward trajec- 
tories were run isentropically from the initialized points for 
14 days. Diabatic descent, which can impact the Antarc- 
tic lower stratosphere over a 14-day period in midwinter, is 
neglected in these calculations. The statistics of the trajec- 
tory analysis are shown in Figure 4. The solid line depicts 
the percentage of trajectories for each day that have passed 
through the ice nucleation criterion (reaching temperatures 
of about 185 K or below). Symbols show the average, min- 
imum and maximum time that the ensemble of trajectories 
spends below the ice frost point (188 K) after the nucleation 
criterion has been met. The temperature statistics suggest 
that the onset of dehydration occurred at the end of June 
when a significant number of the trajectories met the nucle- 
ation criterion for the formation of ice particles. By early 
July, 60 to 90% of the parcels encountered temperatures be- 
low 185 K, suggesting that by this time ice particles nucle- 
ated in nearly all the air parcels initialized w ithin the 188 K 
temperature contour. In addition, the lifetime of ice parti- 
cles nucleated was long enough (5 days or longer) to provide 
sufficient lime for growth and sedimentation of particles to 
lower altitudes. 

5. Model Extinctions 

Next we focus on the extinction measurements from 
CLAES. Previous studies have used these measurements 
to identify PSCs [Mergenthaler et al.^ 1997; Ricaud et a/., 
1995; Massie et ai, 1994]. Here vve use Mie calculations to 
compute aerosol extinction due to ice particles. The calcu- 
lated extinction values are compared to the CLAES observa- 
tions to infer ice particle size and number density. 

Our calculated extinction depend upon temperatures, pres- 
sure, water vapor mixing ratio and fraction of aerosols nucle- 
ated into ice particles. To compute the extinction, we have 
distributed the volume of condensed water (as ice) at equi- 
librium over a fixed and assumed number density of ice par- 
ticles. The size of the ice particles obtained in this manner, 
for a given number density, w^ere used as input parameters 
into a Mie code to convert the volume and size information 
into extinction units [Toon et al., 1990; Massie et ai, 1994). 
The CLAES observations are compared w ith extinction cal- 
culated using different assumptions as to the number density 
of ice particles based on the range of number densities ob- 
served in the in situ Hofmann etal., 1991 1. Maps of the 
calculated extinction are formed by computing the extinc- 


tion for the temperatures at the geographic locations of the 
CLAES measurements when the temperature is below the 
threshold for type II PSCs (188 K). At locations where this 
temperature criterion is not met, the observed CLAES val- 
ues are used in constructing the extinction maps. The maps 
depict what the CLAES instrument would have observed if 
type II PSCs were present in the locations where they are 
predicted to form and if the instrument measured perfectly. 
We examine days from early in the winter since once con- 
densation and ice particle formation begins, the amount of 
water vapor available changes. For the early wintertime 
conditions, we use a mixing ratio of 4.1 ppmv, based on 
the mean values in Figure 1. Figure 5 shows maps of the 
CLAES data and calculated extinction for three different as- 
sumptions of ice particle density. The figure shows that an 
ice particle number density in the range of 10“^- 10“^ cm“^ 
best matches the CLAES observations. The comparison be- 
tween the observed and calculated extinction is good despite 
the uncertainly in the CLAES measurements for high val- 
ues of extinction. Based on the modeling analyses of Jensen 
and Toon [ 1994], the low ice number density implies that the 
clouds were formed in regions with very slow cooling. 

The size of ice particles associated with the concentra- 
tions of 10“^- 10”^ cm““^ is shown in Figure 6. The distri- 
butions are formed from the calculated sizes for the days in 
Figure 5. Ice particle sizes are on the order of 10-30 ^m. 
These should be considered general values since the real at- 
mosphere is likely to contain a distribution of ice particle 
number densities and sizes. In general our inferred ice par- 
ticle number density and size are in good agreement with 
balloon observations over McMurdo station [Hofmann etai, 
1991]. 

6* Comparison of dehydration from POAM 
and MLS 

A recent study of water vapor and extinction data from 
POAM HI examines dehydration in the Antarctic for the SH 
winter of 1998 [Nedoluha et at., 2000]. POAM is a solar oc- 
cultation instrument and makes 14 measurements per day in 
each hemisphere around a latitude circle. The latitude sam- 
pled varies semi-annually as shown in Figure 7. The POAM 
instrument measures water vapor and extinction with a high 
vertical resolution (~1 km) in the lower stratosphere, mak- 
ing it useful for studies of dehydration and PSC occurrence. 
Nedoluha et ai [2000] report dehydration occurring over a 
6 week period starting in mid July with the drop in water 
vapor coincident with the timing of the minimum tempera- 
tures dropping below the ice frost point in air parcels sam- 
pled by the POAM instrument and with the occurrence of 
a significant number of extinction measurements indicating 
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the presence of type II PSCs. 

While the data used in our study cover a different calen- 
dar year, we note that our analysis of the water vapor mix- 
ing ratios from MLS and temperature statistics suggests that 
the onset of dehydration begins prior to early July, and it is 
therefore necessary that we relate our conclusions to those 
of Nedoluha et al. [2000). MLS mixing ratios from 70°S to 
80"^S drop significantly (from 4 to 2.7 ppmv) from mid June 
to July 10 {Figure 1). The area of dehydration and cold tem- 
peratures begins to develop in late June (Figure 3). Vdmel 
et al. [1995] suggest that dehydration and the cold temper- 
ature regions develop initially in the interior of the vortex 
and expand outward. Thus the first observation of dehy- 
drated air depends upon the latitude being sampled. From 
late June through mid July, the POAM instrument observes 
latitudes from 65°S-70°S. The MLS water vapor measure- 
ments in this latitude range do not show evidence of signifi- 
cant dehydration in 10-12 July (Figure 7), whereas substan- 
tial dehydration is evident in the MLS data closer to the pole 
(70°S-80°S) for those dates. Because MLS shows no signs 
of dehydration in mid June and substantial dehydration on 
10 July, we conclude that the onset of severe dehydration 
occurred sometime between mid June and 10 July. Further, 
temperature statistics provide additional support for the on- 
set of dehydration to be around late June. Therefore, if the 
meteorological conditions during the 1 992 and 1998 Antarc- 
tic southern w inters are comparable, then the onset of dehy- 
dration observed by POAM type instruments will always be 
offset as compared to the onset seen by MLS type instru- 
ments. By August, POAM is observing latitudes poleward 
of 70°S and deficits in the water vapor are apparent in the 
data during this period [Nedoluha et al., 2000]. Thus, while 
the analysis of the POAM data provides a much finer ver- 
tical picture and more continuous daily measurements than 
does the MLS data, it is important to note that the latitudinal 
sampling influences the determination of the onset of dehy- 
dration. 

7. Summary 

Examining satellite measurements of water vapor from 
the SH winter of 1992, we find evidence for dehydration due 
to the formation of type II PSCs on potential temperature 
surfaces from 420 to 520 K (16 to 22 km). However, POAM 
observations of water vapor show that dehydration extends 
down to --360 K (12 km). Thus the vertical extent of de- 
hydration is at least 4 km deeper than what MLS can see, 
according to the POAM data. The deficits in water vapor are 
seen in measurements from mid July persisting into Novem- 
ber. Mean mixing ratios decrease by approximately 1 to 2 
ppmv. Temperature statistics suggest that the onset of dehy- 


dration occurs in late June when a large number of cold air 
mass trajectories meet the ice nucleation criteria. Because 
the onset of dehydration initially occurs close to the pole, the 
latitudinal sampling of the atmosphere by different instru- 
ments must be taken into account for the determination of 
the dehydration onset. The development of dehydration over 
the course of the winter shows that initially (June-July obser- 
vations) dehydration is found only in areas which have cold 
temperatures (below 188 K). These areas coincide with those 
with high aerosol extinction. Area calculations show that 
regions which meet the cold temperature criteria and those 
which meet the criteria for dehydrated air expand through- 
out July. By August, the dehydrated regions are larger than 
those of the cold temperatures and high extinction. The de- 
hydrated mixing ratios persist after the elevated extinction 
and cold temperatures diminish. Results from a Mie code 
calculation in conjunction with the extinction measurements 
give a general estimate of aerosol properties. A number den- 
sity between 10““- 10“^ cm“^ best matches the observed ex- 
tinction and, our estimate of ice particle size is on the order 
of 10-30 /xm. 
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